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Abstract. We investigate the non- universal part of the orbital entanglement 
spectrum (OES) of the v — l/i fractional quantum Hall effect (FQH) ground-state 
with Coulomb interactions. The non-universal part of the spectrum is the part 
that is missing in the Laughlin model state OES whose level counting is completely 
determined by its topological order. We find that the OES levels of the Coulomb 
interaction ground-state are organized in a hierarchical structure that mimic the 
excitation-energy structure of the model pseudopotential Hamiltonian which has a 
Laughlin ground state. These structures can be accurately modeled using Jain's 
"composite fermion" quasihole-quasiparticle excitation wavef unctions. To emphasize 
the connection between the entanglement spectrum and the energy spectrum, we also 
consider the thermodynamical OES of the model pseudopotential Hamiltonian at finite 
temperature. The observed good match between the thermodynamical OES and the 
Coulomb OES suggests a relation between the entanglement gap and the true energy 
gap. 
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1. Introduction 

Most condensed matter phases (states of matter) can be characterized using local order 
parameters. However, in some cases, such as systems with topological order, it is not 
possible to obtain a full characterization of a state using only local operators. It has been 
proposed that non-local measurements borrowed from quantum information theory [1], 
such as entanglement, can provide new insights into topological phases phases. The most 
commonly-used quantifier of entanglement, the Von Neumann entanglement entropy, 
measures the entanglement between two blocks of the system. In systems exhibiting the 
Fractional Quantum Hall (FQH) effect, the first topological phases to be experimentally 
realized, such an entanglement measure providing a single number does not provide a 
unique characterization of the many possible states that can occur. A few years ago, it 
was suggested that generalizing the entanglement entropy to the entanglement spectrum 
[2], could reveal much more of the physical properties of a FQH state. 

The notion of entanglement spectrum has now been applied to many different 
systems: Quantum Hall mono-layers [2l|3lll|5l|6l[71|8l[g[T0l[IIl[I3[T3l[Tl[T5], Quantum 
Hall bilayers[l6l [H], quantum spin systems [iHlliniEOlEIlEaElllMlESlEelEg 
[291 l30l EH |32], superconductor |33l l3l] . topological insulators [35l |36l |37] and ultra-cold 
gases |38l |39] . For a system described by a density matrix p, the orbital entanglement 
spectrum (OES) is defined using an angular-momentum orbital decomposition that cuts 
the system into two regions A and B. Such an orbital cut mimics a geometrical cut only 
when the orbitals are localized in space. The reduced density matrix pA is obtained 
by tracing out the B subsystem degrees of freedom, which yields pa = Tr^p. As the 
eigenvalues of pa are non-negative, one can write pa = exp(— H), thus introducing a 
fictitious Hamiltonian "H, whose spectrum is the OES. 

For several FQH states, the count of the low-lying states in the entanglement 
spectrum was numerically shown to match that of the state's edge-theory. Moreover, 
for a realistic Hamiltonian ground state that has a large overlap with a FQH model 
wavefunction, the OES exhibits a low-lying branch with the same entanglement state- 
count as the model wavefunction, but now accompanied by higher entanglement energy- 
levels which were previously not believed to provide useful information about the 
system. In this work, we show that the higher energy levels in the Coulomb OES 
are organized into branches whose structure can be related to virtual particle-hole 
excitations that dress the simpler entanglement spectrum of the model ground-state 
that just characterizes (in its purest form) the topological order of the FQH state. 

This paper is organized as follows. In Section O we introduce the sphere geometry, 
present the concept of the entanglement spectrum for the fractional quantum Hall effect, 
and summarize the main results that have already been obtained in the literature. 
In Section [3], we introduce Jain's composite fermion wavefunctions and their neutral 
excitations. In Section HI we use the composite fermion construction to interpolate 
between Laughlin state and Coulomb interaction ground state entanglement spectrum. 
This interpolation explains the hierarchical structure observed in the OES. In Section 
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O we investigate the behavior of the entanglement spectrum at finite temperature and, 
based on it, conjecture a relation between the energy gap and the entanglement gap. In 
the final section El we present a discussion of these results. 

2. Entanglement spectrum of fractional quantum Hall States on the sphere 

We will consider a system of particles moving on the surface of a sphere (^, cj)) through 
which Nq, magnetic fiux quanta pass. The radius of the sphere is equal to R= ^jN<^/2. 
In the lowest Landau level (LLL), the one-particle orbitals can be expressed as 

1/2 

where (^) is the binomial coefficient, u = cos(6'/2)e"^/^ and v = sin(^/2)e~^'^/^ are the 
spinor variables. These orbitals are eigenstates of Lz, the z-component of the angular 
momentum, with eigenvalues given by / — Nq>/2, where I ranges from to Nq,. The 
orbitals ijji form an approximate real-space partition of the sphere into rings: the north 
pole (resp. south) corresponds to / = Nq> (resp. / = 0). 

Generic fermionic (bosonic) many-body wave functions of N particles and total 
azimuthal angular momentum L*°* can be expressed as linear combinations of Fock 
states in the occupancy basis of the single particle orbitals. Each Fock state can be 
labeled either by A, a partition of into components, or the occupation number 
configuration n{X) = {^^(A), / = A^$, 0}, where n;(A) is the number of times / appears 
in A. We define ^^squeezing" as a two-particle operation on partitions that moves a 
particle from orbital li (resp. I2) to orbital l[ (resp. I2) such that h + h = I'l + I2 
and h < l[ < I2 < h for bosons or li < l[ < I2 < h for fermions. Squeezing defines 
a partial ordering on partitions: if a partition fi can be obtained from a partition A 
using successive squeezing operations, the partition A is said to dominate the partition 
/i (A > yu). Certain model wave functions have a root partition Aq: in their expansion 
on the Fock states, only partitions dominated by Aq can have a non-zero weight. For 
instance, this is the case in the 1/m Laughlin state, in which the occupation number 
configuration of the root partition is given by n(Ao) = {10™'~^10™'~^1 . . .}, where 0™~^ 
denotes m — 1 consecutive empty orbitals. This root partition is ^^{l,m)- admissible" - 
it obeys a generalized Pauli principle which does not allow more than one particle to 
occupy m consecutive orbitals. 

The orbital entanglement spectrum (OES) is obtained by cutting the sphere into 
two parts A and B. Part A contains the Ia first orbitals from the north pole while 
part B contains the Ib = Nq> + 1 — Ia remaining orbitals. Due to the localized nature 
of the Landau level (LL) orbitals, this cut is a reasonable approximation to a spatial 
one [lO]. Tracing over the Ib orbitals gives a reduced density matrix pA with a block- 
diagonal structure. Each block is characterized by two quantum numbers, {Lz,a,Na)- 
the ^-component of the angular momentum for the part-A orbitals and the number of 
particles in part A are the only symmetry-generators present after the partial trace (the 
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For each basis states \^f) E B (i denotes the position in the basis) 
for each basis states I'iff) G A 

create the basis state |^) = |^/) ) 

find k, the index of in the basis describing A® B and the coefficient Ck of 
in the full state |$), including sign from reordering if needed (fermionic case) 
for each index k in the basis of A 

for each index p in the basis of A with p > k 
add to the matrix element p{k,p): Ck * Cp 



Table 1. : Algorithm that computes a given block A of the reduced density matrix 
p = |<i)) ($1 defined by the quantum numbers {L^^a, Na) of the reduced density matrix 
p for the state |$). The complementary space B is defined by {L^^b = L^ — L^^a, = 
N-Na) 



full L symmetry of the original state is lost). For the sphere geometry, the OES is the 
plot of the negative logarithm C,i of pA eigenvalues as a function of Lz^a for a fixed A^^- 
The ^j's are called "entanglement (pseudo)energies" . 

For numerical efficiency we use an algorithm (described in Table [1]) that computes 
each block of the reduced density matrix independently, as for a given pair of quantum 
numbers of the A region, {Lz,a,Na), the only _B-region subspace to consider is that 
defined by {L^^b = — Lz,a,Nb = N — Na)- When we consider a system that is 
not in a pure state, this procedure is repeated for all states present in the density 
matrix. The reduced density matrix thus obtained is then diagonalized using standard 
full diagonalization techniques. 
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Figure 1. Orbital entanglement spectrum for the Laughlin state, (a), and for the 
ground state of the Coulomb interaction, (b) for = 8 fermions, iV$ = 21, Na = 4 
and = 11. A small system-size has been selected for pedagogical purposes. The 
low-lying part of these spectra are almost identical, and exhibit the same structure 
and state-count. In addition to the Laughlin-like branch starting at Lz^a = 24, the 
Coulomb spectrum contains at least two other clearly-defined branches starting at 
Lz,A — 28 and Lz,a — 30. 
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Figure [T] shows a typical OES at z/ = 1/3 for both the Laughhn state (fig. [TJa) 
and the Coulomb interaction ground state (fig. [Hb). In the Laughlin state, the 
entanglement stops at a maximum value of L^^a = L'^'a ■> z component the angular 
momentum of the complementary region B. At this value, a single state is found for 
any values of the total particle number, despite the fact that the corresponding Hilbert 
space dimension of states in A grows exponentially with A^^, the number of particles 
in the sector investigated. When applied to the Laughlin state root configuration, 
the orbital partitioning results in the root configuration 10010010010 for the region A, 
which corresponds to the single state found at the highest possible value of L^^a- While 
the existence of a root partition makes it obvious that no state can be found with a 
higher L^^Ai the fact that there is a unique state for this value is nontrivial and is 
related to a strong constraint on the state decomposition on the Fock basis called the 
"product rule"[lll |17]. The coefficient of a configuration whose two parts after a cut 
can be independently-obtained by squeezing from the state root partition is equal to the 
product of the coefficients of the two disconnected pieces. When decreasing Lz^a from 
this value, the counting of the eigenvalues (the number of eigenvalues at a given L^^a) of 
the model state is the same as the number of levels of the corresponding edge conformal 
field theory (CFT) in the thermodynamic limit; empirically, it consists of counting the 
number of (1, 3)-admissible partitions with the correct angular momentum in I a + AL^ 
orbitals, where AL^ = — Lz^a- This conjecture can be proved using bulk-edge 

correspondence ^2] . In the finite size systems, the equivalence between thermodynamic 
CFT edge counting and orbital entanglement spectrum counting does not hold for all 
Lz^A values. The counting of the spectrum develops finite-size effects that have been 
recently [12] related to the encoding of the Haldane exclusion principle within the model 
state. 

The entanglement spectrum of the Coulomb interaction ground state at filling 
z/ = 1/3 exhibits a branch of low- lying levels displaying the same CFT counting as the 
Laughlin state, separated from higher energy states by an entanglement gap which was 
conjectured to remain finite in the thermodynamic limit [2]. However, this entanglement 
gap closes as Lz^a is reduced, making its definition ambiguous. Using a different Fock 
state normalization in which each LL orbital is normalized by the same factor, it has 
been shown that, in the Coulomb ground state OES, a full gap emerges between low- 
lying states whose state-count is the same as that of the Laughlin state, and higher 
entanglement energy states |1] which had been previously deemed non-universal. In 
this paper (in section H]) we show that even the higher entanglement energy states 
in fact present a universal structure related to particle-hole excitations of the Haldane 
pseudopotential Hamiltonian for the Laughlin state. For now, we notice in the Coulomb 
OES the presence of well separated structures ("christmas-tree"-like branches) with 
states starting at a larger L^^^ than that of the Laughlin-like branch (for example, the 
branches starting at L^'X' = 28 and L^'X' = 30 in Figure [TJb). Understanding these 
branches is the purpose of this paper. 
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3. Composite fermions wavefunctions and their excitations 

Jain's "composite fermion" (CF) picture |13] provides a nice heuristic explanation of 
many features of quantum Hall effect, including the observed incompressible states at 
u = p/{2p + 1) as well as the existence of a compressible the state at z/ = 1/2 (see 
reference for an extensive review on CF). The CF ansatz replaces the strongly 
interacting electrons or bosons by "composite fermions" formed by binding them to n 
flux quanta, where n is even or odd, depending on whether the original "bare" particles 
are bosons {n odd) or fermions (n even). The generic Jain states are given by: 



^'CF = Vlll 



i<j 



(1) 



where Y[i<ji.'^i'^j ~ UjVi)"' binds n flux quanta to each original particle. "Plll is the 
projection operator onto the LLL. is the wave function of the free CF in p effective 
Landau levels, which Jain has called "Lambda levels" (AL). When p such levels are fully 
occupied, equation[T]gives rise to a model state that remarkably accurately approximates 
the ground state at filling u = p/{np + 1). The observed incompressible states at 
v = pI (2p+ 1) can then be thought of as incompressible integer QH CF states u* = p 
where v* is the filling factor of the CFs. Moreover, the compressibility of the v = \/2 
state can be attributed |45| to the formation of a CF Fermi sea state as the effective 
magnetic field felt by the CFs vanishes. In this picture, the u = 1/m Laughlin wave 
function is interpreted as the limiting case n = m — 1 with p = 1 filled CF A levels. A 
schematic view of this picture of the Laughlin state is shown in figure [2^. Except for the 
Laughlin state, the Jain wavefunctions are not known to be unique zero-energy states 
of a model Hamiltonian. However, it has been recently been shown that they do have 
some special "squeezing" properties: the bosonic Jain states at u = p/{p + 1) and the 
bosonic counterpart = f/ ([li<ji'^i'^j ^^e fermionic one at z/ = p/ (2p+ 1) 

have a single root partition and vanish with a power 2 when p + 1 particles are brought 
to the same point [3]. 

The natural way to build quasihole-quasiparticle excitations within the CF 
construction is to assume that the AL's are separated by an effective cyclotron energy 
hu* and to then sort the different excited states with respect to their energy. Thus, the 
lowest energy excited states above the Laughlin ground-state, with energy is hoj*, are 
obtained by exciting one CF in the second AL as shown in figure [2]d. For 2f}Lo* energy 
states, there are two possibilities: two CFs can be put in the second AL (fig. [2t) or one 
CF in the third AL - all other CFs remains in the lowest one (fig. In general, the 
n^^ excited states branch involves the n + 1 lowest A levels. Although this method is 
just a simple phenomenological sketch (as interactions between "composite fermions" 
should surely also be taken into account, and there is no phenomenological reason for 
equal spacing of the levels), the resulting wavefunctions remarkably- well reproduce the 
low energy structure of the Coulomb-interaction Hamiltonian when their variational 
energies are evaluated using the exact Hamiltonian and CF diagonalization|44j . 
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Figure 2. (a) Schematic representation of the Laughhn state as a CF filled lowest A 
level for N = 5 CFs. (b) Creation of a quasiparticle-quasihole excitation above the 
Laughlin state in the ~ sector, (c) and (d) Two different possibilities to create 
an excited state with energy 2huj*. 



Monte Carlo methods have been used extensively to quantitatively explore the 
predictions of the CF picture[46], proving useful in computing quantities such as 
predicted energies of CF states and their overlaps with states obtained by exact 
diagonalization. However, calculations of entanglement spectra requires a high accuracy 
in the decomposition of the wave function on the Fock basis. With the exception of 
several special cases [3], Monte Carlo techniques fail to reach the accuracy needed and 
we needed to implement a new exact projection method. The first step of this method 
is to use the Jack recursion formulas of [HI HT] to expand Y[j<ki'^j'^k — VjUk)"' (which 
is a Jack polynomial) in terms of symmetric monomials or slater determinants (for n 
even or odd respectively). The product with is computed explicitly, using only the 
resulting wave function symmetry properties to reduce the computation time. Finally, 
the projection is applied and only terms fully in the LLL are kept. We found that 
the clustering properties characterizing how "composite particle" quasihole-quasiparticle 
excitation states vanish as particles come together are determined by the number of the 
A level involved: the bosonic states and the bosonic counterpart of the fermionic states 
(fermionic states divided by a vandermonde determinant) constructed from the /c'th 
AL vanish when k + 1 particles are brought to the same point as the second power of 
the difference in the particle coordinates. Therefore the excited states also exhibit a 
non-trivial root configuration. In all the cases we have studied, we noticed the following 
properties of states of energies less or equal to 2hw*: such excited states can be realized 
in a number of different ways: all CFs in the lowest AL, one or two CFs in the second 
AL or one CF in the third AL. The last possibility does not introduce any new states 
compared to the ones that are produced using the first two ways. To obtain this result, 
the LLL projection is crucial as it creates linear dependencies between the projected 
states and greatly reduces the number of independent states. For instance, for = 6 
fermions and = 15 in the L*°* = sector, there are 51 states with 2 CFs or less in 
the second AL and none in the third LL whereas there are 57 states whose energy is less 
than or equal to 2hL0*; after projection and re-orthonormalization these two numbers 
are reduced to 36 and the states spaces spanned by these two sets of states are identical. 
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4. Prom Laughlin to Coulomb entanglement spectrum 

In this section, we use the admixture of virtual quasiparticle-quasihole excited states into 
the ground state to interpolate between the Laughlin state and the Coulomb interaction 
ground state at the same filling factor, and examine the effect of this on the orbital 
entanglement spectrum. As both states have L = we only consider the zero angular 
momentum sector . The z/ = 1/3 Laughlin state is known to be the densest zero-energy 
eigenstate of the Haldane pseudopotentials[lH] by Vi = 1 and Vn>i = 0. Using this 
interaction, we could obtain the low energy L = states and use them to reconstruct 
the Coulomb interaction ground state. However, these states do not have an exact model 
wavefunction construction that would lead to a simple understanding of the additional 
structure in the orbital entanglement spectrum of the Coulomb state. In particular 
the L = excited states of this model interaction do not feature a non trivial root 
configuration, contrary to the CF approach. Therefore, we first use the CF quasihole- 
quasiparticle excitations states whose construction was described in the previous section 
to iteratively construct a basis of L = states (the iteration being the number of AL 
levels occupied) whose effective cyclotron energy is less than a given effective cyclotron 
energy. We approximate the Coulomb ground state using: 

'.nhui* 



CF I coulomb j 



CF 



(2) 



- CF 



where the sum runs over all the L = states of the iterative basis which spans the 
space of states whose effective cyclotron energy is less than nhu*. This state is then 
normalized, which translates into a global shift of the entanglement spectrum. As the 
effective cyclotron energy is increased, the partition that dominates all the basis states 
changes and successively dominates the previous basis (see Table [2]). The L™^^''s of 
the successive root partitions correspond to those of the additional structures observed 
on Coulomb ground state entanglement spectrum - thereby allowing us to predict the 
starting point of the higher entanglement energy branches. In the Table [21 the number 
of L = states, the topmost partitions and the corresponding L^j^ values for the 
accessible effective cyclotron energy are given for = 8 and = 21. 



n 


number of L = states 


Topmost partition 


r max 


1 


1 


10010010010|01001001001 


24 


2 


4 


11000100100|00100100011 


28 


3 


8 


11000100100|01000010011 


28 


4 


14 


11001001000|00010010011 


30 



Table 2. Characteristics of the iterative basis which spans the space of states whose 
effective cyclotron energy is less than nhui* for iV = 8 and iV$ — 21. The last column 
indicates the value corresponding to the topmost partition with respect to the 

cut defined by La — 11 and Na — A - shown by a vertical line in the partitions. The 
total Hilbert space has 31 L = states - we are using at most 14 of them. 
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1(a) n = 2 
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(b) n = 3 



Figure 3. Orbital entanglement spectrum of different \^n)- n = 2 (a), n = 3 (b) 
and n = 4 (c) for iV = 8 fermions, Ng, = 21, Na = 4 and Ia — H- The structures 
observed on the Coulomb state OES progressively appear as n is increased and we 
obtained completely seperatcd branches when n is even. The last observed branch in 
the Coulomb spectrum, for which L™^^ = 32, is expected to be described by an higher 
effective cyclotron energy CF wavefunctions. 



The comparison of the orbital entanglement spectrum for the first three interesting 
|\E'„) displayed in Figure [3l and the Coulomb state spectrum (Fig. [1]) shows that the 
CF quasihole-quasiparticle excitation basis allows a step-by-step reconstruction of the 
different Coulomb entanglement spectrum structures. Given the method used, it is 
obvious that the spectra of |\E'„) and l"^ coulomb) would get closer as the effective cyclotron 
energy is increased; in that case, we include more and more AL's and, in a finite- 
size calculation, we start diagonalizing in larger and larger parts of the full Hilbert 
space. However, nontrivially, the structures of the Coulomb ground state entanglement 
spectrum are sorted by the effective cyclotron energy, in the sense that the first, second, 
etc. branches of the Coulomb spectrum can be obtained by only using CF states with 
n = 1,2, etc. We checked these properties for the fermionic u = 1/3 Laughlin state 
for up to = 8 particles and for the bosonic u = 1/2 one for up to = 10 particles 
(Fig. H]). The different scales in the eigenvalues of the Coulomb interaction ground 
state density matrix seem to be linked to the different effective cyclotron energy scales 
involved. 




Figure 4. Orbital entanglement spectrum of for the Coulomb interaction ground state 
(a) and different |*„): n = 2 (b) and n = 4 (c) for TV = 10 bosons, 7V$ = 18, Na ^ 5 
and Ia = 9. 



Hierarchical structure in the orbital entanglement spectrum in Fractional Quantum Hall systemslO 



5. Orbital entanglement spectrum at finite temperature 

We now try to reproduce the Coulomb ground state entanglement spectrum using a 
thermal density matrix. For a system at finite temperature T, the density matrix is 
given by: 



P 



iexp(-/3if) 



(3) 



where /3 = 1/T, Z = Tt [exp(— and H denotes the Hamiltonian of the system. In 
our case, we want to compare the entanglement spectrum of the Coulomb-interaction 
ground state to the entanglement spectrum of the density matrix describing finite- 
temperature corrections to the Laughlin state and hence the simplest choice is to take 
H to be the pseudopotential interaction Hamiltonian. At T = we recover the pure 
Laughlin-state entanglement results. In contrast to the previous section, we now have to 
take into account not only the L = states, but states with all possible values of L, and 
all the Lz sectors of the pseudopotential Hamiltonian. This involves a very large number 
of states and Equation [3] can be realized exactly only for small systems. By choosing 
an ad-hoc "entanglement temperature" temperature T, the entanglement spectrum of p 
and the one of the Coulomb ground state can be made very similar as shown on figure 
[5l Although the full shape of the spectrum is very similar, we notice that the Coulomb 
spectrum exhibits degeneracies which are not present in the thermal density matrix 
approach. For example, the rightmost levels are degenerate for the Coulomb-interaction 
ground state and are no longer degenerate when the thermal density matrix is used, 
despite its rotational-invariance. It can be analytically shown that these degeneracies 
are linked to the L = states that enter the density matrix. When the sum of Equation 
[3] is restricted to just the L = states, these degeneracies are recovered (fig. [5b). 



;a) Coulomb 



10 12 



(b) |3=7, all L 




Figure 5. Orbital entanglement spectrum for the Coulomb ground state (a), for p 
defined by equation [3] with /3 = 7 (b) and p' defined by equation [3] but with the sum 
restricted to the L = states with /? = 3 (c) for N = 6, = 15, = 3 and La = 8. 
The value of the inverse temperature /3 has been chosen so that the spectra are similar 
to that of the Coulomb ground state. This naturally leads to a choice of two different 
temperature values for p and p' as the number of states considered in these two cases 
is very different. The general shape of these spectra is very similar. 



To analyze bigger systems, we consider a simpler model. In the thermal density 
matrix, we only kept the two lowest energy branches of the incompressible state: the 
ground state wave function and the magneto-roton mode. This is the single-mode 
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approximation of the entanglement spectrum. The dispersion relation of the magneto- 
roton mode is that obtained through by the Vi interaction. The magneto-roton states 
I^L^f )' obtained by putting one CF in the second AL, consist of A^ — 1 multiplets whose 
total angular momentum ranges from L = 2 to L = N. Even though a L = 1 state is 
naively expected, this multiplet is systematically suppressed by the LLL projection, as 
pointed out in |19]. Thus we use the following approximation for p: 



P 



(4) 



where 



LI. 



L is the total angular momentum of the multiplet states and is their energy with 
respect to the pseudopotential interaction Hamiltonian. We then compute the orbital 
entanglement spectrum associated with this density matrix for various /3 values and for 
up to = 10 fermions (see Fig. E]). 
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Figure 6. Entanglement spectrum of the single-mode approximation p given by 
equation [4| for a system of = 8 fermions and ~ 21 flux quanta, (a): j3 = 0, 
(b): /3 = 1, (c): (3 = 1, (d): (3 = 12. The cut used is iV^ = 4 and Ia = H. 



In every case we have studied, when the temperature is infinite, the universal 
Laughlin-like part of the entanglement spectrum cannot be distinguished. As 
temperature is decreased, two branches of levels split: a lower one whose state-count is 
the same as that of the Laughlin state, and an upper one whose average "entanglement 
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(see figure E]), to define an entanglement gap [2], given by Sl,^^ = ^n+i.L,,^ - ^n,L,,A 
where ^n,L^ a highest entanglement energy at L^^a that belongs to the Laughhn-state 
structure. For example, in the infinite temperature limit in figure El this can be done 
for Lz^A = 24,23,22, at which the Laughlin-like levels of counting 1,1,2 respectively 
are separated from the higher energy states by a small but visible gap. Moreover, the 
upper-branch state-count, and the Lz^a value at which it starts, are the same as those 
of the first branch after the Laughlin one in the Coulomb ground state entanglement 
spectrum. As the magneto-roton mode can be obtained from excitations to the first AL 
level, only the first entanglement spectrum above the Laughlin one can be fitted in this 
single-mode approximation. 

To more precisely characterize the behavior of the two branches as a function of the 
temperature, we calculated the entanglement gap 6 function of /3 = 1/T. The 

entanglement gap, shown in Figure [71 decreases as the temperature is increased from 
T = but, starting from a certain temperature value of the order of the energy gap, it 
reaches a plateau. This behavior is the same for all the cases we studied, independent 
of the number of fermions. Moreover, to see how our crude single-mode approximation 
affects the entanglement spectrum, we considered the next energy-level feature, taking 
all CF excitations states whose effective cyclotron energies are less than or equal to 
2tvjj*. These spectra are presented in Figure [H where we observe the emergence 
of a third branch in the thermal density matrix, matching the third branch of the 
Coulomb entanglement spectrum. The main behavior as the temperature is varied is 
not substantially affected by these additional states. However, it should be noticed that 
an additional structure appears and that the entanglement gap at high temperature no 
longer exhibits a plateau but just a change of slope (shown in the Figure [7]). 
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Figure 7. Entanglement gap J^max for different temperature values and different 
number of particles: iV = 8 in red, TV = 9 in blue and = 10 in green and in 
turquoise = 8 with all the states whose effective cyclotron is less than 2fiu!*. The 
vertical line is the temperature corresponding to the true energy gap. 
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Figure 8. : Entanglement spectrum of p given by a formula completely analogue to|4] 
in which the 2fia;* effective cyclotron energy states are taken into account, for a system 
of TV = 8 fermions and iV* 21 flux quanta, (a): /3 = and (b): /3 = 10. The cut 
used is Na — 4 and Ia = H- 



6. Discussion 



The entanglement spectrum can be used to determine the universahty class of realistic 
Hamiltonians in a topologically ordered phase: by identifying the state-count of the 
lowest-lying entanglement branch of the spectrum with the state-count (characters) of 
the CFT of an edge theory, we can in principle predict that the ground-state of the 
system lies in a certain topological phase. In this paper we showed that, for the case of 
the z/ = 1/3 Coulomb interaction, not only the low-lying branch (with the Laughlin-state 
state-count), but also the higher' "entanglement-energy" branches exhibit a nontrivial 
characteristic structure, related to dressing of the simple model ground state that gives 
rise to the lowest branch of the entanglement spectrum by zero-point fluctuations of 
particle-hole collective excitations, when the difference between the model Hamiltonian 
and the Coulomb Hamiltonian is added back as a "perturbation" . We have explicitly 
showed that the higher-energy branches correspond to the wavefunctions for adding 
particle-hole excitations on top of the Laughlin ground-state. The correspondence 
exhibits quantum number matching {L^^a and counting of the levels in a specific branch) 
with the Coulomb spectrum if the model wavefunctions we use for the excitations 
are Jain's CF wavefunctions. We then performed a single-mode approximation of the 
Coulomb spectrum by calculating the reduced thermal density matrix of the Laughlin 
state augmented by the magnetoroton mode. We found that by this method we 
could obtain the first branch above the Laughlin-state branch in the entanglement 
spectrum of the Coulomb state. This exercise also appears to supports the idea that 
the entanglement spectrum of the ground-state of a realistic Hamiltonian contains 
information not only about the universality class of the ground-state but also about 
its excitations, which in a generic Hamiltonian (unlike in free fermion Hamiltonians 
FQH model Hamiltonians) will be represented in the ground state properties though 
zero-point fluctuations of collective modes. 
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